The Mycoplasma conjunctivae genome sequencing, annotation and analysis

BACKGROUND: The mollicute Mycoplasma conjunctivae is the etiological agent leading to infectious keratoconjunctivitis (IKC) in domestic sheep and wild caprinae. Although this pathogen is relatively benign for domestic animals treated by antibiotics, it can lead wild animals to blindness and death. This is a major cause of death in the protected species in the Alps (e.g., Capra ibex, Rupicapra rupicapra). METHODS: The genome was sequenced using a combined technique of GS-FLX (454) and Sanger sequencing, and annotated by an automatic pipeline that we designed using several tools interconnected via PERL scripts. The resulting annotations are stored in a MySQL database. RESULTS: The annotated sequence is deposited in the EMBL database (FM864216) and uploaded into the mollicutes database MolliGen http://cbi.labri.fr/outils/molligen/ allowing for comparative genomics. CONCLUSION: We show that our automatic pipeline allows for annotating a complete mycoplasma genome and present several examples of analysis in search for biological targets (e.g., pathogenic proteins).

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

Bioinformatics tools and databases for the study of human growth hormone

Advances in novel molecular biological diagnostic methods are changing the way of diagnosis and study of metabolic disorders like growth hormone deficiency. Faster sequencing and genotyping methods require strong bioinformatics tools to make sense of the vast amount of data generated by modern laboratories. Advances in genome sequencing and computational power to analyze the whole genome sequences will guide the diagnostics of future. In this chapter, an overview of some basic bioinformatics resources that are needed to study metabolic disorders are reviewed and some examples of bioinformatics analysis of human growth hormone gene, protein and structure are provided.

Molecular cytogenetics and gene mapping in sheep (Ovis aries, 2n = 54)

The development of a completely annotated sheep genome sequence is a key need for understanding the phylogenetic relationships and genetic diversity among the many different sheep breeds worldwide and for identifying genes controlling economically and physiologically important traits. The ovine genome sequence assembly will be crucial for developing optimized breeding programs based on highly productive, healthy sheep phenotypes that are adapted to modern breeding and production conditions. Scientists and breeders around the globe have been contributing to this goal by generating genomic and cDNA libraries, performing genome-wide and trait-associated analyses of polymorphism, expression analysis, genome sequencing, and by developing virtual and physical comparative maps. The International Sheep Genomics Consortium (ISGC), an informal network of sheep genomics researchers, is playing a major role in coordinating many of these activities. In addition to serving as an essential tool for monitoring chromosome abnormalities in specific sheep populations, ovine molecular cytogenetics provides physical anchors which link and order genome regions, such as sequence contigs, genes and polymorphic DNA markers to ovine chromosomes. Likewise, molecular cytogenetics can contribute to the process of defining evolutionary breakpoints between related species. The selective expansion of the sheep cytogenetic map, using loci to connect maps and identify chromosome bands, can substantially contribute to improving the quality of the annotated sheep genome sequence and will also accelerate its assembly. Furthermore, identifying major morphological chromosome anomalies and micro-rearrangements, such as gene duplications or deletions, that might occur between different sheep breeds and other Ovis species will also be important to understand the diversity of sheep chromosome structure and its implications for cross-breeding. To date, 566 loci have been assigned to specific chromosome regions in sheep and the new cytogenetic map is presented as part of this review. This review will also summarize the current cytogenomic status of the sheep genome, describe current activities in the sheep cytogenomics research sector, and will discuss the cytogenomics data in context with other major sheep genomics projects.

A Staphylococcus xylosus isolate with a new mecC allotype

Recently, a novel variant of mecA known as mecC (mecA(LGA251)) was identified in Staphylococcus aureus isolates from both humans and animals. In this study, we identified a Staphylococcus xylosus isolate that harbors a new allotype of the mecC gene, mecC1. Whole-genome sequencing revealed that mecC1 forms part of a class E mec complex (mecI-mecR1-mecC1-blaZ) located at the orfX locus as part of a likely staphylococcal cassette chromosome mec element (SCCmec) remnant, which also contains a number of other genes present on the type XI SCCmec.

Novel pseudo-staphylococcal cassette chromosome mec element (ψSCCmec57395) in methicillin-resistant Staphylococcus pseudintermedius CC45

Genetic characterization of methicillin-resistant Staphylococcus pseudintermedius (MRSP) from Thailand and Israel revealed the presence of a predominant atypical clonal lineage which was not typeable by SmaI-PFGE and SCCmec typing. All the atypical isolates (n = 34) belonged to CC45 (30 ST45 and 2 ST179 isolates, 1 ST57 isolate, and 1 ST85 isolate). The isolates originated from healthy and diseased dogs and cats, as well as from the environment of one clinic. Cfr9I-pulsed-field gel electrophoresis (Cfr9I-PFGE) and dru typing permitted the further distinction of CC45 isolates from the two different countries. Microarray analysis identified genes that confer resistance to β-lactams (mecA; blaZ), aminoglycosides [aac(6')-Ie-aph(2')-Ia; aph(3')-III; ant(6)-Ia], macrolides and lincosamides [erm(B)], tetracyclines [tet(M)], trimethoprim [dfr(G)], streptothricin (sat4), and chloramphenicol (catpC221). Fluoroquinolone resistance was attributed to specific amino acid substitutions, i.e., Ser84Leu in GyrA and Ser80Ile and Asp84Asn in GrlA. A novel pseudo-staphylococcal cassette chromosome (ΨSCCmec57395) element was identified in MRSP strain 57395 (sequence type ST45) by whole-genome sequencing. The 12,282-bp ΨSCCmec57395 element contained a class C1 mec gene complex but no ccr genes. In addition to the methicillin resistance gene mecA, ΨSCCmec57395 also carried determinants of resistance to heavy metals, such as arsenic, cadmium, and copper. Bsu36I restriction analysis of the ΨSCCmec57395 element amplified by long-range PCR revealed the presence of ΨSCCmec57395 in the 33 additional isolates of MRSP CC45. The ΨSCCmec57395 element represents a new class of SCCmec and has been identified in MRSP of CC45, which is a predominant clonal lineage in Israel and Thailand.

A point mutation in cpsE renders Streptococcus pneumoniae nonencapsulated and enhances its growth, adherence and competence

BackgroundThe polysaccharide capsule is a major virulence factor of the important human pathogen Streptococcus pneumoniae. However, S. pneumoniae strains lacking capsule do occur.ResultsHere, we report a nasopharyngeal isolate of Streptococcus pneumoniae composed of a mixture of two phenotypes; one encapsulated (serotype 18C) and the other nonencapsulated, determined by serotyping, electron microscopy and fluorescence isothiocyanate dextran exclusion assay.By whole genome sequencing, we demonstrated that the phenotypes differ by a single nucleotide base pair in capsular gene cpsE (C to G change at gene position 1135) predicted to result in amino acid change from arginine to glycine at position 379, located in the cytoplasmic, enzymatically active, region of this transmembrane protein. This SNP is responsible for loss of capsule production as the phenotype is transferred with the capsule operon. The nonencapsulated variant is superior in growth in vitro and is also 117-fold more adherent to and more invasive into Detroit 562 human epithelial cells than the encapsulated variant.Expression of six competence pathway genes and one competence-associated gene was 11 to 34-fold higher in the nonencapsulated variant than the encapsulated and transformation frequency was 3.7-fold greater.ConclusionsWe identified a new single point mutation in capsule gene cpsE of a clinical S. pneumoniae serotype 18C isolate sufficient to cause loss of capsule expression resulting in the co-existence of the encapsulated and nonencapsulated phenotype. The mutation caused phenotypic changes in growth, adherence to epithelial cells and transformability. Mutation in capsule gene cpsE may be a way for S. pneumoniae to lose its capsule and increase its colonization potential.

Global phylogenomic analysis of nonencapsulated Streptococcus pneumoniae reveals a deep-branching classic lineage that is distinct from multiple sporadic lineages.

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated Streptococcus pneumoniae (Non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if Non-Ec-Sp evolve sporadically, if they have high antibiotic non-susceptiblity rates and a unique, specific gene content. Here, whole genome sequencing of 131 Non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multi-locus sequences types ST344 (n=39) and ST448 (n=40). All ST344 and nine ST448 isolates had high non-susceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic Non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic Non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic Non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P=0.005). In contrast, sporadic Non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, Non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of pneumococcal conjugate vaccines, Non-Ec-Sp may become more prevalent.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Looking the cow in the eye: deletion in the NID1 gene is associated with recessive inherited cataract in Romagnola cattle

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.

The novel macrolide-Lincosamide-Streptogramin B resistance gene erm(44) is associated with a prophage in Staphylococcus xylosus.

A novel erythromycin ribosome methylase gene, erm(44), that confers resistance to macrolide, lincosamide, and streptogramin B (MLSB) antibiotics was identified by whole-genome sequencing of the chromosome of Staphylococcus xylosus isolated from bovine mastitis milk. The erm(44) gene is preceded by a regulatory sequence that encodes two leader peptides responsible for the inducible expression of the methylase gene, as demonstrated by cloning in Staphylococcus aureus. The erm(44) gene is located on a 53-kb putative prophage designated ΦJW4341-pro. The 56 predicted open reading frames of ΦJW4341-pro are structurally organized into the five functional modules found in members of the family Siphoviridae. ΦJW4341-pro is site-specifically integrated into the S. xylosus chromosome, where it is flanked by two perfect 19-bp direct repeats, and exhibits the ability to circularize. The presence of erm(44) in three additional S. xylosus strains suggests that this putative prophage has the potential to disseminate MLSB resistance.

Bacillus Calmette-Guérin strain differences have an impact on clinical outcome in bladder cancer immunotherapy

BACKGROUND Whether the commonly used bacillus Calmette-Guérin (BCG) strains Connaught and Tice confer different treatment responses in non-muscle-invasive bladder cancer (NMIBC) is unknown. OBJECTIVES To compare clinical efficacy, immunogenicity, and genetics of BCG Connaught and Tice. DESIGN, SETTING, AND PARTICIPANTS A prospective randomized single-institution trial with treatment of 142 high-risk NMIBC patients with BCG Connaught or Tice. INTERVENTION Patients were randomized to receive six instillations of BCG Connaught or Tice. For experimental studies, BCG strains were compared in C57Bl/6 mice. Bladders and lymphoid tissues were analyzed by cytometry and the latter cultivated to detect live BCG. BCG genomic DNA was sequenced and compared with reference genomes. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Recurrence-free survival was the primary end point of the clinical study. The Kaplan-Meier estimator was used for estimating survival and time-to-event end points. Nonparametric tests served for the analysis of the in vivo results. RESULTS AND LIMITATIONS Treatment with BCG Connaught conferred significantly greater 5-yr recurrence-free survival compared with treatment with BCG Tice (p=0.0108). Comparable numbers of patients experienced BCG therapy-related side effects in each treatment group (p=0.09). In mice, BCG Connaught induced stronger T-helper cell 1-biased responses, greater priming of BCG-specific CD8(+) T cells, and more robust T-cell recruitment to the bladder than BCG Tice. Genome sequencing of the BCG strains revealed candidate genes potentially involved in the differential clinical responses. CONCLUSIONS BCG strain may have an impact on treatment outcome in NMIBC immunotherapy. PATIENT SUMMARY We compared the efficacy of two commonly used bacillus Calmette-Guérin (BCG) strains for the treatment of NMIBC and found that treatment with BCG Connaught prevented recurrences more efficiently than BCG Tice. Comparison of the immunogenicity of the two strains in mice indicated superior immunogenicity of BCG Connaught. We also identified genetic differences that may explain the differential efficacy of the Connaught and Tice BCG strains. TRIAL REGISTRATION NCT00003779.